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Abstract. We present the angular correlation function of the X-ray population of 1063 XMM- 

CO [ Newton observations at high Galactic latitudes, comprising up to ^30000 sources over a sky area 

of ~125 sq. degrees in the energy bands: soft (0.5-2 keV) and hard (2-10 keV). This is the largest 

sample of serendipitous X-ray sources ever used for clustering analysis purposes to date and the 

results have been determined with unprecedented accuracy. We detect significant clustering signals 

\^ . in the soft and hard bands (^10c7 and ^5a, respectively). We deproject the angular- coiTelation 

(-H ' function via Limber's equation and calculate the typical spatial lengths. We infer that AGN at 

— ' redshifts ^1 are embedded in dark matter halos with typical masses of logM ^ 12.6/i^' Mq and 

lifetimes in the range ^3 — 5 x 10^ years, which indicates that AGN activity is a transient phase in 

O ' the life of galaxies. 
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1. INTRODUCTION 

^ . Active galactic nuclei (AGN) are the brightest persistent extragalactic sources known, 

O I with their X-ray emission the most common feature among them. Thanks to their 

IT) ■ large bolometric output, AGN can be detected through cosmological distances, which 

makes them essential tracers of galaxy formation and evolution, as well as the large- 
scale structure of the Universe. Clustering studies of AGN at redshift ~1, when strong 
structure formation processes took place, are key tools for understanding the underlying 
mass distribution and evolution of cosmic structures. 
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2. THE X-RAY DATA 

In this work we use the sample presented in [9]. The selected XMM-Newton/EPIC-pn 
observations fulfilled the following criteria: 1. High galactic latitude fields i\b\> 20°) in 
order to minimize the contamination from Galactic sources; 2. Fields with at least 5 ks 
of clean exposure time; and 3. Fields free of bright and/or extended X-ray sources. 

If there were observations carried out at the same sky position, the overlapping area 
from the observation with the shortest clean exposure time was removed. The resulting 
sample comprised 1129 observations. For the purposes of this work we have also re- 
moved the observations belonging to the Virgo Cluster, M31, M33, Large Magellanic 
Cloud and Small Magellanic Cloud fields, ending up with a final sample of 1063 obser- 



vations. The overall sky coverage of the sample is 125.52 deg^ comprising 31288 and 
9188 sources in the soft (0.5-2 keV) and hard (2-10 keV) bands, respectively. 

3. THE ANGULAR CORRELATION FUNCTION 

The two-point angular correlation function w{6) determines the joint probability of 
finding two objects in two small angular regions 5Q.i and 5^2 separated by an angular 
distance with respect to that of a random distribution ([11]). 

To calculate the angular correlation function we have used the estimator proposed 
by [7] w{ei) = DD-2DR+RR ^ ^fJ\lQXQ DD, DR and RR are the normalised number of 
pairs of sources in the z-th angular bin for the Data-Data, Data-Random and Random- 
Random samples, respectively. To produce the random source sample against which we 
have compared the real source sample searching for overdensities at different angular 
distances, we have tried to mimic as closely as possible the real distribution of the 
detection sensitivity of the survey. The method employed here is extensively described 
in [1] and [4]. The errors in different angular bins are not independent from one another. 
To estimate the errors we have followed the co variance matrix method described in [10]. 

The angular correlation function can be described by a power-law model in the form 

(\ 1— y 
^ j , where 1 — 7 is the slope and Oq is the angular correlation length. 

We have fitted the data using a x^ technique. 

4. RESULTS 

In the soft band (0.5-2 keV) we detect a high-significance (~10a) clustering signal 
with a correlation length of Oq = 22.9 ±2.0 arcsec and a slope of 7— 1 = 1.12 ±0.04 
after correcting for the integral constraint. If we ignore that correction in our fits, the 
correlation lengths are comparable within the error bars but the power-law becomes 
significantly steeper (7— 1 = 1.29 ±0.04). 

In the hard band (2-10 keV) the power-law becomes steeper (7— 1 = 1.33) and 
the clustering is stronger (Oq = 19.T^^'q arcsec) although marginally consistent with 
the results in the soft band within the icr error bars. The clustering detection is still 
very significant (~5(7). The sources detected in the hard band are less biased against 
absorption, and if the unified model of AGN is correct (the obscuration of the nucleus 
is due to orientation effects only) one might not expect significant differences in the 
clustering properties of obscured and unobscured sources. However, accurate angular 
clustering measurements in this band have been difficult because of the limitations 
caused by the small-number statistics (e.g. [5], [1], [13]). 

5. INVERSION OF LIMBER'S EQUATION 

The two-dimensional angular correlation function is a projection in the sky of the real 
three-dimensional spatial correlation function t, (r) along the line of sight, where r is 
the physical separation between sources (typically in units of h^^ Mpc). The angular 





FIGURE 1. Angular correlation function in the 0.5-2 keV (left panel) and 2-10 keV (right panel) bands. 
Solid dots are the observed data. Triangles represent the averaged w{9) for a random dataset. Overplotted 
is the best-fit x^ with and without fixed slope (dashed and solid lines, respectively). 



amplitude Gq can be related to the spatial amplitude ro by inverting Limber's integral 
equation ([12]). For this, we assumed a redshift selection function for our sample 
derived from the best-fit luminosity-dependent density evolution (LDDE) model of the 
X-ray luminosity function of [3] for all bands under study. Hence, we found values of 
ro =12.25 ±0.12 and 9.9 ± 2.4 h^^ Mpc in the soft and hard bands, respectively, for 
a clustering model constant in comoving coordinates, while for clustering constant in 
physical coordinates we obtained ro =6.54 ±0.06 and 5.7 ±1.4 h^^ Mpc, respectively. 

6. DARK MATTER HALOES AND THE LIFETIME OF AGN 



The spatial clustering values derived above can be used to estimate the mass of the dark 
matter haloes (DMH) in which these sources are embedded. A commonly used quantity 

for such an analysis is the bias parameter, that is usually defined as b^{z) = f^fg). 

where ^agn{^iZ) and ^dmh{^,z) are the spatial correlation functions of AGN and DMH 
evaluated at 8 h^^ Mpc, respectively. The former value has been calculated in this 
work whereas the latter can be estimaded using the expressions in [11]. In a A-CDM 
cosmology, the bias is function of both the mass of DMH and redshift ([8]). From 
this we derive an average mass of (logMfyMH) — 12.60 ±0.34 h^^ Mq. This result 
is fully in agreement with that of [6] who found logMoMH = 12 A— 12.8 h^^ VIq in 
the COSMOS survey, although it is worth noticing that our results come from a mostly 
unidentified X-ray sample. Similarly, if we compare our results with the clustering study 
of [2] in the AEGIS field we see that our mean DMH mass is intermediate between that 
of blue and red galaxies, which is in agreement with the location of X-ray selected 
AGN in the color-magnitude diagram. We can estimate the lifetime of AGN using 
the mean DMH masses calculated above and making some simple assumptions. We 
have followed the method proposed by [8], assuming that we are sampling the most 
massive DMH at a given redshift z and that each DMH hosts an active AGN at any 



given time. Hence we have tAGN{z) = tu{z) 
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where ?(/(z) is the Hubble time at 



redshift z, 4>(z) is the comoving density of AGN for which we used the predictions of the 



luminosity function of [3], and ^dmh{z) is the comoving density of DMH that can be 
estimated following the Press-Schechter approximation. At z ~ 1 and for haloes larger 
than logMmin = 12.6 h^^ M©, we get (i>DMH — 2 x 10"^ h^ Mpc"-^. This yielded to an 
AGN duty cycle in the range tACN/tu = 0.054 — 0.078. In our cosmological framework, 
the Hubble time atz= 1 is ~5.8 Gyr. The estimated lifetime of AGN is hence in the 
range Iagn = 3.1 — 4.5 x 10^ yr. 

7. CONCLUSIONS 

We have studied the angular correlation function of a large sample of serendipitous X- 
ray sources from 1063 XMM-Newton observations at high Galactic latitudes in the 0.5- 
2 keV (soft), and 2-10 keV (hard) energy bands, respectively. Our sample comprises 
31288 and 9188 sources in the soft and hard bands, respectively, covering ~125.5 deg^ 
in the sky, thus being the largest sample ever used in clustering investigations. 

We found significant positive angular clustering signal in the soft (~10(7) and hard 
(~5(7) bands. The result in the hard band clears up the debate on whether X-ray sources 
detected in this band cluster or not, since a number of past works had reported different 
inconclusive results ranging from a few o detections to no detection at all. 

We inverted Limber's equation, assuming a given redshift distribution for our sources, 
in order to estimate typical spatial correlation lengths. We used these values to calculate 
the rms fluctuations of the AGN distributions within a sphere of radius 8 h^^ Mpc, and 
compared them with that of the underlying mass distribution from the linear theory in 
order to estimate the bias parameter of our X-ray sources. The bias depends on the mass 
of the dark matter haloes (DMH) that host the AGN population. From the computed bias 
values we have estimated a typical DMH mass of {logMoMn) — 12.60 ±0.34 h^^ Mq. 

The typical AGN lifetime derived from the Press-Schechter approximation at redshift 
z ~ 1 lies in the range in the range tACN = 3.1— 4.5 x 10^ yr. This interval is significantly 
shorter than the time span between that redshift and the present thus suggesting the 
existence of many AGN generations, and that a significant fraction of galaxies may 
switch from a quiescent phase to AGN activity, and vice versa, several times throughout 
their lifes. 
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